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We study the scattering properties of topological solitons on obstructions in the form of holes and 
barriers. We use the 'new baby Skyrme' model in (2+1) dimensions and we model the obstructions 
by making the coefficient of the baby skyrme model potential - position dependent. We hnd that that 
f--^ ' the barrier leads to the repulsion of the solitons (for low velocities) or their complete transmission 

(at higher velocities) with the process being essentially elastic. The hole case is different; for small 
velocities the solitons are trapped while at higher velocities they are transmitted with a loss of 
energy. We present some comments explaining the observed behaviour. 
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PACS numbers: 11.10.Lm,12.39.Dc,03.75.Lm 



> ■ I. INTRODUCTION 

m . 

Consider a moving particle encountering a potential barrier or a potential hole. Then, the barrier slows down the 
particle and if the energy of the particle is too low relative to the barrier's height, the particle is reflected. In the 
' potential hole case the particle speeds up in the hole and is always transmitted. 

However, consider now the same process in quantum mechanics; then in both cases we have reflection and trans- 
mission, the relative magnitudes of which depend on the parameters of the potential and on the energy of the particle. 
Recently, some work has been done on the scattering properties of solitons of the nonlinear Schrodinger model 
was pointed out in Ref. |5j that the scattering of solitons at low velocities resembles a classical particle 
in the sense that the soliton maintains its integrity and follows a well-defined trajectory, with the difference that, 
nevertheless, the soliton can be reflected by a potential hole. Such reflection is not possible in classical mechanics. 

Motivated by these results, we have decided to look at this problem in the case of topological solitons asking ourselves 
the question as to whether this case will resemble more the classical particle case or the quantum mechanical systems? 

To look at this more carefully we have chosen to study this problem on the example of a (2+1) dimensional system 
so that we could see the effects of the transverse direction and to allow for the radiation waves, if they are generated, 
to escape more easily. The model we have chosen is the baby Skyrme model discussed in detail in many publications 
0, Q • The mutual interactions of Skyrmions in scattering events have been a topic of long standing interest and are 
summarized in Refs. 0,0,^^- Experimental realizations of topological solitons are possible in spinor Bosc-Einstein 
condensates [llHl^ . 

The Lagrangian density of the "baby Skyrme" model contains three terms, from left to right: the pure S 2 sigma 
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model, the Skyrme and the potential terms: 



c = d„<t> ■ d»4> - 9s [(d^ ■ d^) 2 - (d^4> ■ cW)(c"V • d v (t>)\ - V{4>) (1) 

where 

V(ct>)=fi(l-cj>l). (2) 
The vector </> lies on the unit sphere S 2 hence <fi ■ <fr = 1. 

To have a finite potential energy the field at spatial infinity is required to go to 03 = ±1, 0i = 02 = 0. In this work 
we choose "the vacuum" to be defined as 03 = +1. 

Note that our boundary condition has defined a one-point compactification of R2 , allowing us to consider on the 
extended plane i?2 U 00 topologically equivalent to S 2 . In consequence, the field configurations are maps 

$:S 2 — >S 2 . (3) 

which can be labelled by an integer valued topological index Q: 

Q = 7~ J {p x $ x ®v£) dxd V- (4) 

As a result of this non-trivial mapping the model has topologically nontrivial solutions which describe "extended 
structures" , which have been called baby skyrmions. A soliton is then the simplest field configuration corresponding 
to Q = 1 which minimises the total energy and which can be calculated from C by taking all terms in it with 
positive signs. This soliton field configuration has to be found numerically; however, as discussed in previous papers 
this problem reduces to having to solve an ordinary differential equation for a "profile function" /(r), where r 
is the radial distance from the position of the soliton. 

Note that the soliton is exponentially localised and its asymptotic behaviour is controlled by /i. 

In the next section we describe a possible way of introducing an "obstruction" into our model and discuss the results 
of our studies. In Sec. Ill we present the results of our numerical simulations. Sec. IV presents some concluding 
remarks. 



II. POTENTIAL OBSTRUCTION 

There are various ways of introducing a potential hole or a potential barrier. However, given that the soliton field, 
strictly speaking, is never zero, even though it vanishes exponentially as we move away from its position, this potential 
has to be introduced in such a way that it does not change the "tail" of the soliton i. e. it has to vanish when 03 = 1 . 
A possible way to do this is to add an extra term to the Lagrangian which vanishes when 03 = 1. Of course, there 
are many possible choices of such terms but given that our Lagrangian already contains a term with such a property 
we exploit this fact and choose to add a(l — <j>\) in some region of x and y. We choose this term to be independent 
of y so that the obstruction on the potential energy landscape, located in some finite region of x, say at positive x, 
resembles a trough in the "hole" case or a dam in the "barrier" case. Then sending the soliton from a point well away 
from this obstruction, i.e. initially placed at some sufficiently negative x, in the positive x direction, we can study 
the effects of the obstruction. 

In our numerical simulations we have chosen the obstruction to be constant in a small range of x; this effectively 
corresponds to taking /x in the original Lagrangian to be given by /iq for x in the range of the obstruction and y,\ 
elsewhere. 
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The case when fi > Mi corresponds to a barrier (dam), and when /z < Mi we have a hole (trough). We have 
performed many numerical simulations of such systems, varying both the sign and value of /io — fJ-i and the velocity of 
the incoming soliton. We have also checked that, initially the solitons were far enough from the obstruction so that 
the incoming solitons can be considered to be free. 

In the next section we present some results of our simulations. Both the cases of the hole and of the barrier have 
been studied and, as we shall see, they have produced very different results. The scattering by the barrier was found 
to be elastic, with hardly any radiation being generated in the process; by comparison the motion through the trough 
was inelastic with an interesting pattern of the decrease of velocity of the transmitted soliton. 

III. NUMERICAL SIMULATIONS 
A. General comments 

We have performed most of our simulations on a 350x350 square grid with lattice spacing being 0.15. Thus the 
lattice extended from -26.5 to +26.5 in each direction. The soliton was initially placed at x = —12, y = 0. Its size 
was determined by the choice of parameters, 8 S and fi which were chosen as 9 S = 0.5 and [i\ — 0.6. This produced 
a soliton which was essentially localised to about 40 lattice points in each direction i.e. to a region of 6 x 6 in real 
space. Thus, placing the soliton at x = — 12 there was no problem with any boundary effects. 

Next we put our obstruction from x = to x = 30 * 0.15 = 4.5. We considered various values of the height of the 
obstruction. As they all led to qualitatively similar results we performed most of our simulations for — /i = 0.5 
and —0.1. The same is true when we varied the width of the barrier although changing the width changed the value 
of v cr . 

All our simulations were performed using a 4th order Runge-Kutta method with the field 4> being rescaled every few 
iterations. The time step of our simulations was taken to be 0.001. We used fixed boundary conditions and later we 
used also the absorbing boundary conditions at the edges of the lattice. This we generated by successively decreasing 
the magnitude of at the last 5 rows and columns of the lattice. 

To generate the time dependence we calculated the initial configurations with the soliton located at x = —12 and 
x = —11.9 and we defined as being proportional to the difference between these two 4>' orthogonalised with respect 
to (f> at x — —12. Varying the constant of proportionality we changed the initial velocity. This is an approximate way 
of introducing the initial condition corresponding to a moving soliton, which is a very good approximation for small 
values of velocities, when the relativistic effects are negligible. Using this method we have to calibrate the velocity 
- to check directly what the initial velocity is. Another way to proceed would involve introducing the correct time 
dependence into the initial ansatz (with Lorcntz factors etc) and then calculating ^ directly from this expression. 
Given the fact that the initial configurations are calculated numerically and involve some extrapolations we chose 
the easier option mentioned above. Our results show that our procedure was exceptionally good; the moving soliton 
varied its height very little and it did not generate any perceptible amount of radiation. 

B. Barrier 

First we have considered the effects of a barrier. Hence, between x = and x = 4.5, wc put \i\ — /io = —0.1. Then 
we placed the initial soliton at various values of x and calculated its energy without altering the shape of the soliton. 
This has told us what the barrier is like as seen by the soliton. The calculated energy plot is shown in Fig. 1 

We see that, although the original barrier is in the shape of a square barrier, the soliton perceives the barrier as 
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FIG. 1: Energies of our basic soliton as a function of its position 



a smooth hump or hill. This is, of course, due to the finite size of the soliton. As seen in Fig. 1, the barrier effects 
stretch from x ~ — 4 to x ~ 10. 

The soliton placed at x = —12 should be far away from the effects of the barrier. Next we performed a series 
of simulations placing the soliton at x = —12 and sending it towards the barrier at various speeds. We have found 
extremely elastic behaviour in all cases. For low velocities (v < v cr ~ 0.2427) the soliton bounced off the barrier 
while for v > v cr it was transmitted. At velocities close to v cr the soliton slowly climbed the ascending slope of the 
effective barrier of Fig. 1, either getting over it (v > v cr ) or falling back. In all the cases the process was elastic. This 
was seen through the plots of energy density, where we did not see any significant radiation energy. Also the values 
of the final velocity of the soliton were essentially the same as the initial ones, indicating that no significant energy 
loss had taken place. To calculate the initial velocity we used the values of the time when the soliton reached x = — 5 
on the way towards the barrier. To calculate the final velocity we used the times it reached x = —5 and x — — 11 for 
the reflected soliton and x = 11 and x = 21 for the transmitted one. This way we were calculating the velocities in 
the region where there were no effects of the barrier. In Fig. 2 we plot the modulus of the outgoing versus incoming 
velocity. 

We note that the curve is essentially a straight line suggesting a simple linear relationship. In fact, the values of 
the incoming and outgoing velocities differ from each other in 3 or 4th decimal points. So we can treat them as equal, 
within the numerical accuracy of our procedure. 

The value of the critical velocity can be estimated by observing the trajectory of the soliton; if the soliton makes 
it to the 'top of the barrier', which is around x ~ 2.25, then it goes through, otherwise, it is reflected. In Fig. 3 we 
present plots of the time dependence of the soliton position in two cases; one corresponding to the velocity of just 
below the critical value, namely, v ~ 0.2426 and one just above at v ~ 0.2428. 




FIG. 3: Time dependence of the position of the soliton initially placed at x = 



— 12: left v < v cr , right v > v, 
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FIG. 4: Energies of our basic soliton as a function of its position 



We note that, in each case, the slopes of the curves before and after the soliton interacted with the barrier, are 
the same. Thus we see that the scattering is completely elastic; essentially no energy is lost during the climb of the 
barrier and because of this the velocity is unchanged, to the degree of accuracy of our calculation. 

Looking at the static configuration with /i = 0.7 we note that its energy is 1.968 which is exactly the energy that the 
incoming soliton has to have to be able to transfer all its kinetic energy into the potential energy - to be able to climb 
the barrier and to be transmitted. Thus v cr is determined by the energy that the incoming soliton corresponding to 
fi = 0.6 has to have so that its energy corresponds to the energy of the static soliton of \i = 0.7. 

Note that one can give the following nonrelativistic argument which supports our claims. First note that the mass 
of the soliton is given by the total energy (c = 1) so M ~ 1.912. Then ^f- ~ 1.912 * 0.05890329/2 = 0.05631 which 
is in agreement with the difference of energies SE — 1.968 — 1.1912 ~ 0.056. 

This assumes an elastic behaviour and this is what we have seen in our numerical simulations. Thus in the case of 
the bump - the whole kinetic energy of the incoming soliton is converted into the potential energy of the soliton at 
the top of the barrier and then the soliton can be transmitted elastically. 



C. Hole 



We placed the hole in the same place as the barrier and we took /xi — /io = 0.5. Next, as in the barrier case, we 
placed the initial soliton at various values of x and calculated its energy. This has told us what the hole is like as 
seen by the soliton. The calculated energies are shown in Fig. 4 Once again we see that the effective hole as seen by 
the soliton is quite smooth. We have performed many simulations and have found that, again, there exist a critical 
velocity v cr . Above v cr the soliton is transmitted and below it falls in and becomes trapped in the hole. 
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FIG. 5: Time dependence of the position of the soliton initially placed at x = —12: left v < v cr , right v > v cr 

The critical velocity is around v cr ~ 0.155. Solitons started off from x = —12 fall into the hole and stay there 
oscillating and gradually losing their energy by emitting radiation. Two typical trajectories are shown in Fig. 5. The 
left picture shows a trapped soliton with initial velocity v ~ 0.09, ie below v cri and the right one shows a transmitted 
soliton with initial velocity v ~ 0.16. 

For velocities above v cr we have a transmission, but this time the outgoing soliton is much slower. In Fig. 6 we 
present the plot of the outgoing velocities as a function of the incoming ones. 

We note not only the significant decrease of the velocity but also the oscillations in the values of the outgoing 
velocities. 

We have looked at the details of the scattering but we have not succeeded in revealing the origin of the oscillations. 
We have rerun some simulations with absorption of the waves at the boundaries; the results were essentially the same 
suggesting that the origin of the oscillations has nothing to do with any waves of emitted radiation bouncing off the 
boundaries. 

Thus - the origin of the effect is associated with the interaction of the soliton with the hole and the emission of the 
radiation during this process. To check this we have looked at the time spent in the hole and its relation to the time 
needed for the soliton to traverse the hole had the soliton moved with the initial or the final velocity. 

In Fig. 7a we plot the time spent by the soliton in the hole, as a function of the incoming velocity. In Fig. 7b we 
plot three curves; all as a function of the initial velocity. They present the difference of the time the soliton spent in 
the hole from which is subtracted the time needed to traverse 1 the hole with the initial velocity, the same with the 
final velocity and the same with the average velocity. 

We see a dramatic difference in all three expressions: 

It is clear that the time in the hole is related to the oscillations of the velocities. It is also clear that as the incoming 
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FIG. 6: Velocities of transmitted solitons, x axis - incoming, y axis - outgoing 

velocity is much larger than the outgoing velocity the expression for the time spent in the hole - the time needed to 
traverse the hole at a given velocity is the lowest when we take the largest velocity, i.e., the incoming velocity. It is 
also clear that the potential hole increases the velocity - so all our expressions in Fig. 7b are negative. 

Looking at the curves in Fig. 7b we see that the oscillations in outgoing velocity are pretty much reproduced by the 
time spent by the soliton in the hole and almost eliminated when we take the middle curve in Fig. 7b. This suggests 
that during the motion through the hole the soliton velocity changes in a complicated way (gradually producing the 
oscillations). If, in the plots in Fig. 7b we use the initial velocity the oscillations show up in the difference of times; if 
we use the final velocity we overcompensate and we have the opposite oscillations; taking the average velocity accounts 
for most of the oscillations leaving a curve which is much smoother. This curve, however, does show the extra velocity 
in the hole, which of course is related to the incoming velocity. 

Looking at the curve we note that its features are as expected; we note a slow rise with with the increase of Vi nc , 
which is consistent with the loss of energy in the hole. 

Note that this time the process is very different from the situation of a barrier. In the barrier case the incoming 
soliton has some kinetic energy which it has to convert to the potential one to rise up the barrier. This seems to be 
done clastically, with essentially no radiation. In the hole case the soliton has a kinetic energy and as it moves into 
the hole it gains extra energy as well. This time the extra energy is partially converted into the increased speed of 
motion but some of it is got rid of via radiation. In fact, this is what we have seen in our numerical simulations which 
did show some radiation. Then, when it reaches the other side of the hole and it starts 'climbing out of the hole' the 
soliton may not have enough energy to get out. Hence, for low speeds it gets trapped in the hole. 

Of course, we do not really understand why the soliton does not convert all its extra energy into the increased 
kinetic energy and why the velocities give us the very interesting patterns that can be seen from the plot in Fig. 6. 
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FIG. 7: left - time spent in the hole; right - time in the hole minus the time needed to traverse the hole with various velocities. 
Top curve - initial velocity, middle average, bottom - outgoing one 

Clearly there must be an internal mode of the soliton that gets excited and which radiates the excess energy. 
This point requires further study. 

IV. CONCLUSIONS 

We have looked at a system involving a topological soliton in two dimensions and a potential, of both a barrier and 
a hole-type. 

When the soliton was sent towards the barrier its behaviour resembles that of a particle. Thus at low energies the 
soliton was reflected by the barrier and at higher energy it was transmitted. The scattering process was very elastic. 
During the scattering the kinetic energy of the soliton was gradually converted into the energy needed to 'climb the 
barrier'. If the soliton had enough energy to get to the 'top' of the barrier then it was transmitted, otherwise it slid 
back regaining its kinetic energy. 

Note that the soliton size is related to the parameters of the model and so depends on [i. Hence, during the climb 
of the barrier, the soliton altered its size (it decreased a little) - to fit the local value of /i; when it got through, or 
slipped back, its size returned to it original value. This is what one would expect in an elastic scattering and this 
is what we saw in the numerical simulations. In fact, the soliton size oscillated a little, around its 'correct' value 
and the amplitude of these oscillations has not changed much during the scattering process and the final oscillations 
resembled the original ones. 

In the hole case, the situation was very different. This time, the soliton gained an extra energy as it entered the 
hole. Some of this energy was converted into kinetic energy of the soliton, some was radiated away. So when the 
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FIG. 8: The density of kinetic energy in a simulation just after the soliton has passed the hole. T = 150, Vi n = 1.9165 

soliton tried to 'get out' of the hole it had less kinetic energy than at its entry and, when this energy was too low it 
remained trapped in the hole. During the scattering process, like in the case of a barrier, the soliton size changed and 
its oscillations increased significantly. Afterwards they stayed like this - with much higher amplitude of oscillations 
than before. Hence the increase in oscillations is related to the inelasticity of the process and the emitted radiation. 

The scattering of a topological soliton on a hole is thus reminiscent of a classical particle under the influence of 
friction. We have not found any indication of nonclassical reflection as in the case of Nonlinear Schrodinger solitons 
However, in contrast to a classical particle with simple velocity-dependent friction, we have found interesting 
oscillations hinting at some underlying resonant mechanism. 

We have looked in detail at the behaviour of the solitons during their scattering process and, so far, have not found 
a satisfactory explanation of the observed 'oscillations' in the outgoing velocities (in the case of transmission by the 
hole). 

Clearly, this is related to the properties of the radiation - the only plausible explanation we can find is that the 
radiation is sent out in separate bursts. This is very much what we saw in the actual simulations. In Fig. 8 we present 
a plot of the density of the kinetic energy seen in a typical simulation. 

Looking at the plot we see distinctive waves of energy moving along the 'trough' (i.e. in the y directions). The 
'high' peak is the kinetic energy density of the moving soliton. 
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Clearly we need to find an explanation of the observed radiation pattern. But, this lies outside the scope of this 
paper. 
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